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In this work we investigate the inclusive and exclusive photoproduction of dileptons, which is 
relevant for the physics programme to be studied in the proposed electron-proton collider, the 
LHeC. In the inclusive case, the process is sensitive to the parton distribution functions in the 
photon whereas the exclusive channel is connected to the small- £ QCD dynamics. For the latter, 
we investigate the role played by saturation physics at a very high energy scenario. The estimates 
for production cross sections and the number of events are presented. 

PACS numbers: 13.60. Fz,13.90.+i, 12. 40.-y,13.60.-r,13.15.Qk 



INTRODUCTION 



The parton distributions functions (PDFs) of the pho- 
ton are fairly determined by current experiments. It is 
an open subject the issue of generating them from avail- 
able high energy data and evolving them using perturba- 
tive QCD approach. It is timely to investigate their im- 
plications for inclusive deep inelastic processes induced 
by photonic parton distributions which can be tested by 
the forthcoming high energy beams experiments. One of 
them is the Deep Inelastic Electron-Nucleon Scattering 
at the LHC (LHeC) machine, which is an extension of the 
LHC experiment as an ep-collider [l| . This proposed ex- 
tension will open a new kinematic window - the photon- 
proton centre-of-mass energy can reach up to TeV scale, 
which it is a very nice region for small- a; physics. Being 
specific, the energy of the incoming proton is delivered 
by the LHC beam, and a list of possible scenarios are 
considered for the energy of incoming electron as E p = 7 
TeV and E e — 50 — 200 GeV, corresponding to the center 
of mass energies of y/s = 2y/E^E~ e ~ 1.18 - 2.37 TeV 0. 
The planned integrated integrated luminosity is of order 
10 — 10 2 fb _1 . Here, we shall consider inclusive dilepton 
(Drell-Yan) production in high energy ep collisions which 
are mostly sensitive the the quark content of photon. It 
is a complementary process to the heavy-quark produc- 
tion in ep reactions that is sensitive to the gluon content 
of photon. 

On the other hand, the exclusive photoproduction of 
lepton pairs, 7 TV -> t*(-> £ + £~) + N (with N = p, A) 
can be studied using well established high energy factor- 
ization approaches. As an example, using simple models 
for the elementary dipole-hadron scattering amplitude 
that captures the main features of the dependence on 
atomic number A, on energy and on momentum trans- 
fer t the differential cross section for exclusive dilepton 
photoproduction can be estimated [3j. Such an inves- 
tigation is complementary to conventional partonic de- 
scription of timelike Compton scattering, which consid- 
ers quark handbag diagrams at leading order in a s and 
simple models of the relevant generalized parton distri- 



butions (GPDs). These calculations are as well input 
in electromagnetic interactions in pp and AA collisions 
to be measured at the LHC. They can also be studied 
with beams of protons or anti-protons as recently done 
by the Tevatron machine. For instance, the CDF Collab- 
oration studied the exclusive production of muon pairs, 
pp — » pp + at lower invariant masses Q. 

The aim of this work is twofold - first, we show pre- 
dictions for the photoproduction of dileptons at planned 
LHeC energies within high energy factorization schemes. 
The inclusive photoproduction cross section including the 
resolved and direct contributions is calculated, focusing 
mostly on the dilepton invariant mass distribution. We 
then move to the dilepton exclusive photoproduction, an- 
alyzing the production cross section using the color dipole 
picture which give us a simple picture of process at small- 
x regime. The present work is complementary to the cal- 
culation presented in Ref . [j| , where the photoproduction 
of massive gauge bosons was investigated as a tool to test 
the Standard Model and at same time to gain insight on a 
very specific physics beyond SM, namely the WW-f cou- 
plings. There, the production of W , Z was computed 
and predictions were made for the number of events pro- 
duced at the LHeC regime. This article is organized as 
follows. The basic formulae to calculate the inclusive and 
exclusive (diffractive) photoproduction of virtual photons 
are presented in next section. Our numerical results for 
cross section and event rates are presented in section Hill 
followed by the correspondent discussion. The summary 
and and conclusions are presented in section HVl 



II. INCLUSIVE AND EXCLUSIVE DILEPTON 
PHOTOPRODUCTION 



Let start by considering the Drell-Yan inclusive cross 
section at next-to-leading order for the process 7p — > 
£ + £~ + X. In particular, we will use the NLO matrix el- 
ements and parton distribution functions (PDFs) in the 
A/5 scheme. The differential cross section on the dilep- 



2 



ton mass is given by [6j, 
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where one has r = M 2 /s. The quantities x 7 and x p 
are the parton momentum fractions from the photon and 
proton, respectively. Here, q 1 , g 7 are the quark and gluon 
distribution functions for the photon whereas q p , g p are 
the corresponding functions for the proton. The first 
term in equation above contains the dominating leading 
order quark-antiquark initial state and the associated real 
and virtual corrections. The second term includes the 
g p + q 1 and q p + g 1 initial states whereas the last term 
is the Compton term. In the numerical calculations we 
set the factorization and renormalization scales at the 
invariant mass of the dilepton pair, /i 2 = Q 2 = M 2 . The 
quantities Pi(x p , a; 7 , Q 2 ) are written down as, 
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The calculation above can be compared to the updated 
calculation for the massive gauge boson photoproduction 
done in Ref. [f|. They include the direct and resolved 
photon contributions to the production rates in an anal- 
ysis designed to the LHeC machine. 

As the inclusive dilepton photoproduction is con- 
cerned, a better experimental signal can be obtained con- 
sidering their exclusive production. This is particularly 
true if one takes high photon luminosity as typical for 
coherent interactions in heavy ion colliders or e + e _ col- 
liders. The high luminosity and the very large centre-of- 
mass energy to be available by the LHeC machine open 
a window for large rates of exclusive dileptons in that 
case. In the high energy limit it is related to the time-like 
Compton scattering, where a virtual photon is produced 
in the photoproduction process. Such a process is better 
viewed in the target rest frame, commonly named the 
color dipole approach. Let us now to introduce the rele- 
vant equations in the color dipole picture for the timclike 
Compton scattering. In this formalism 0] , the scattering 
process jp — > j*p is assumed to proceed in three stages: 
first the incoming real photon fluctuates into a quark- 
antiquark pair, then the qq pair scatters elastically on 



the proton, and finally the qq pair recombines to form 
a virtual photon (which subsequently decays into lepton 
pairs). The amplitude for production of the exclusive 
virtual photon final state is given by [1, 

A^' p (x,Q,A) = V V / d 2 r [ dz** h - h (r,z,Q) 
f h.h J Jo 
x A qq (x,r,A)t> h - h {r,z,0), (5) 

where ^^(r, z, Q) denotes the amplitude for a photon 
to fluctuate into a quark-antiquark dipole with hclicitics 
h and h and flavor /. The quantity A q q(x,r, A) is the 
elementary amplitude for the scattering of a dipole of size 
r on the proton, A denotes the transverse momentum lost 
by the outgoing proton (with t = —A 2 ), x is the Bjorken 
variable and Q 2 is the photon virtuality. Here, we are 
considering the photoproduction case and then only the 
transversely polarized overlap function contributes. For 
a given quark flavour /, when summed over the quark 
helicities it gives, 
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where we have defined the quantities e\ 2 = 
and z = (1 — z). Accordingly, the photon virtualities are 
Q\ = (incoming real photon) and Q 2 = —Q 2 (outgoing 
virtual photon). For simplicity, we consider space-like 
kinematics. In Ref. [lCj, the complete time-like calcu- 
lation is performed in the k± -formalism and the results 
are somewhat larger than the presented here. Therefore, 
the current estimation is a lower bound for the exclusive 
dilepton protoproduction. 

Taking the imaginary part of amplitude, which is dom- 
inant at high energies, the elastic diffractive cross section 
is then given by, 
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In our numerical calculations we consider saturation 
models which correctly describe exclusive processes at 
high energies like vector meson production, diffractive 
DIS and DVCS. Our baseline model will be non-forward 
saturation model of Ref. Q . Its advantage is to give the 
t dependence of elastic differential cross section without 
the necessity of considerations about the impact parame- 
ter details of the process and the overall normalization is 
determined without assumptions about the elastic slope. 
The non-forward scattering amplitude can be written as: 

A qq -{x, r, A) = 2irR 2 p e- B|t| N (rQ sat , x) , (8) 

with the asymptotic behaviors Q 2 at (x, A) ~ 
max(QQ, A 2 ) exp[— Aln(x)]. Specifically, the t de- 
pendence of the saturation scale is parametrised as 
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in order to interpolate smoothly between the small and 
intermediate transfer regions. The form factor F(A) = 
exp(— B\t\) catches the transfered momentum depen- 
dence of the proton vertex, which is factorised from the 
projectile vertices and does not spoil the geometric scal- 
ing properties. Finally, the scaling function N is obtained 
from the forward saturation model [Tl| , whose functional 
form resembles the analytical solution for the nonlin- 
ear QCD equations in the asymptotic energy regime. 
Here, we investigate the exclusive photoproduction of a 
heavy timelike photon which decays into a lepton pair, 
7P — > £ + £~p. Therefore, for the £ + t~ invariant mass 
distribution from the virtual 7* decay we have (with 
Q 2 = M 2 +e _), 
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In the next section we compute the numerical results 
for the inclusive and exclusive photoproduction of dilep- 
tons focusing on LHeC regime of energy /luminosity. We 
also investigate the sensitivity to to the theoretical un- 
certainties for both processes. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

Our starting point is a study of the inclusive dilepton 
photoproduction for the proposed LHeC machine [2, 0] ■ 
Using the design with a electron beam having labora- 
tory energy of E e = 70 GcV, the center of mass energy 
will reach E cm — W lp = 1.4 TeV and a nominal lumi- 
nosity of order 10 33 cm _2 s _1 . In Fig. [1] is presented 
the differential cross section , Eq. (TTJ), as a function of 
the dilepton invariant mass (with Mg+i- > 3 GeV/c). 
Here, we have summed the resolved and direct contri- 
butions and through the calculations proton structure 
functions of CTEQ6M [l2[ and photon structure func- 
tions of GRV [13j have been used with p 2 — M 2 +i _ . 
For sake of comparison it is also shown the result for 
the DESY-HERA energy range. When comparing the 
distinct energies, deviations are larger for higher invari- 
ant masses. The behavior is similar for small invariant 
mass and the overall normalization differs by a factor 
about ten. Our estimates for the inclusive DY photo- 
production cross sections are the following. One gets 
Omc(7 + P ^ £ + £~X) ~ 0.78 nb for the integrated cross 
section with M e + e - > 3 GeV/c. This results can be com- 
pared to the HERA energy, which gives an integrated 
cross section of 0.12 nb. These are rough estimates once 
we have not introduced the theoretical uncertainty com- 
ing from using distinct PDFs for proton and photons. 
The uncertainty concerning the variation of hard scale is 
small in a NLO QCD computation. 

We also compute the electron-proton total cross sec- 
tions for mass cut M e + e - > 3 GeV/c and the correspond- 
ing number of events. The latter has been computed us- 
ing N ev = a(ep — > £ + £~ + X)Li nt . At this point we 



FIG. 1: Differential cross sections for the inclusive photopro- 
duction of dileptons as a function of dilepton invariant mass 
for DESY-HERA and LHeC energies. 



consider the acceptance in the dilepton channel as 100%. 
The photoproduction cross section is calculated by con- 
voluting the Weizsacker- Williams spectrum 
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with the differential hadronic cross section. One has 
Qmin = mly/il-y) and we impose a cut oiQ 2 min = 0.01. 
An integrated luminosity L int of 10 fb _1 0] is assumed 
in order to compute the number of events, N ev . The final 
result is a(ep -> + X) = 11 pb for the LHeC en- 

ergy regime and the number of events is around 1.1 x 10 5 . 
The number of events is considerable as expected for an 
inclusive process. 

Let us now move to the exclusive dilepton photopro- 
duction. In the following we compute the integrated 
cross section, performing the integration over \t\ < 1 
GeV 2 . It should be noticed that the |i|-dependence is 
distinct from other implementations of the dipole cross 
section. However, for integrated cross sections, these de- 
viations do not play an important role. In Fig. [5] we 
focus on the invariant mass (M 2 +£ „ > 1 GeV 2 /c 2 ) de- 
pendence for fixed values of energy, as computed from 
Eqs. ([7]) and (jTOJ) . Once again we show the results for 
the LHeC energy and DESY-HERA for comparison. Ac- 
cordingly, the spectrum is dominated by small invariant 
masses of order M 2 +e _ < 20 GeV 2 , independent of en- 
ergy. In the saturation models, the qualitative behav- 
ior can be obtained considering dominance of small size 
dipoles configuration. The forward amplitude reads as 
A oc (Q 2 at /M 2 ) times a logarithmic enhancement on 
M 2 . Therefore, qualitatively the differential cross sec- 
tion behaves like da/dM 2 oc (M 2 )- S [l + log(M 2 )], with 
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S m 3. The total cross section, (j(jp — >• £ + £~ p), inte- 
grated over dilepton invariant mass Mg+£- > 1.5 GeV 
for energies W » 10 2 GeV can be adjusted by a power 
fit cr = 3pb(W/W ) 0A6 , with W = 1 GeV. The inte- 
grated cross section gives cr eX c(7 + P — > £ + £~ + p) — 80 
pb with the cut M e + e - > 1.5 GeV/c and also one has 
<j{ep — > £ + £~ + p) = 2.1 pb for the LHeC energy regime 
and the number of events is around 2 x 10 4 . The number 
of events is one order of magnitude smaller than for the 
inclusive case. This is consistent with the HERA trend 
of a 10% ratio for diffractive over inclusive processes. 
The extrapolation of current study for the electropro- 
duction case is straightfoward as within the color dipole 
approach the deeply virtual Compton cross section is in 
general following the geometric scaling property. That is, 
acompton(£, Q 2 ) = CTcompton (t) , where the scaling vari- 
able Q is given in terms of incoming photon virtuality 
and the corresponding saturation scale, r = Q 2 / Q 2 aat {x) . 

It is timely to discuss here the limitation of the space- 
like kinematics considered in current calculation. In such 
a case, the wavef unction for an incoming photon with 
spacelike virtuality q 2 — —Q 2 < is very well known 
and can easily used to estimate the exclusive diffractive 
cross section, taking Q 2 = M 2 +e _ for the virtuality of 
final state photon. Within the color dipole approach, the 
correct timelike kinematics has been addressed in Ref. 
fTj ] in the scope of exclusive diffractive photoproduction 
of heavy gauge bosons, like Z° production. However, 
a phenomenological study for timelike Compton scatter- 
ing has not been addressed yet. The main difficulty is 
the numerical computation of scattering amplitude as 
the wavefunction for a outgoing photon with timelike 
q 2 — M 2 +£ _ > leads to an overlap function which gives 
an integrand strongly oscilatory as a function of trans- 
verse dipole size r as noticed in [141 ]. On the other hand, 
the calculation of timelike Compton scattering using the 
fc_i_-factorization formalism has been done in Ref [l(| as 
discussed in previous section. It was found the the ratio 
of the cross section i?(timelike/spacelike) is q 2 -dependent 
and is of order 3 or 4 at high energies. 

We notice that the non-forward amplitude is assumed 
to be partially factorized concerning the x and t depen- 
dences, Eqs. (8) and (9), as the saturation scale also de- 
pends on the momentum transfer t. The model we are us- 
ing describes correctly the ^-dependence of vector meson 
electroproduction and deeply virtual Compton scattering 
(DVCS) as well. In general, within the dipole approach 
the extension of amplitude to finite momentum transfers 
is set to be completely factorized. This assumption, due 
to its simplicity, is widely used in the literature, although 
it contradicts the Regge-pole theory, as was shown e.g. in 
studies of DVCS or general parton distributions (GPD), 
as shown e.g. in Refs. [15j, |l6|. We believe that the dis- 
tinct assumptions for the i-dependence will introduce an 
uncertainty on the overall normalization as we are con- 
sidering only ^-integrated cross sections in current study. 
On the other hand, this issue would be important when 
distributions on momentum transfer is concerned. 
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FIG. 2: (Color online) Differential cross sections for the ex- 
clusive photoproduction of dileptons as a function of dilepton 
invariant mass for DESY-HERA and LHeC energies. 



Finally, let us compare the present calculation to pre- 
vious studies. In Ref. @ the inclusive dilepton photopro- 
duction has first addressed focusing on the DESY-HERA 
energy regime. The current calculation does an upgrade 
for those results using the state of art concerning new 
parton distribution functions for the proton and photon. 
Our predictions for the LHeC are consistent with an en- 
ergy extrapolation from the HERA machine. However, 
the LHeC machine is probing the PDFs at smaller val- 
ues of parton momentum fraction, x p , than for HERA. 
For instance, for invariant mass of 4 GeV/c, one has the 
minimum momentum fraction probed i p = 8x 10~ 6 for 
LHeC against i p = 4x 10~ 4 for HERA. Similar trend 
occurs for the photon parton function as x p — > 1 and 
x 1 goes to M 2 1 s. Therefore, LHeC will be an outstand- 
ing place to study the small- x physics in addition to the 
PDFs for photons. Concerning the exclusive dilepton 
photoproduction, our estimates can be compared to re- 
cent studies in literature. For instance, in Ref. [l(| the 
exclusive dilepton photoproduction is computed within 
the fc_L-factorization formalism. We have checked that 
the invariant mass distribution is consistent with ours 
with deviations only on the overall normalization. In 
particular, the extrapolation for the LHeC energy gives 
an integrated cross section around 1 nb (using the cut 
M 2 +£ _ > 1.5 GeV 2 /c 2 ). This is very large compared to 
our estimates and the origin of deviation is the fact that 
the correct time-like configuration produces a larger cross 
section compared to the space-like configuration we have 
considered. This issue is discussed in detail in Ref. 
and the main point is that the time-like calculation is 
a factor 3 or 4 compared to the space-like photon con- 
figuration. In addition, the timelike Compton scattering 
was also investigated in Ref. [17| within the general- 
ized parton distributions (GPD) factorization approach. 
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The latter is a robust study as the interference at the 
amplitude level with the pure QED subprocess is also 
computed. For the energies we are considering here the 
QED Bcthe-Heitler contribution saturates around 28.5 
pb for the integrated cross section whereas the timelike 
process contributes with 15-25 pb. Disregarding interfer- 
ence contribution (in high energies they are comparable 
to timelike and QED contributions), the total cross sec- 
tion is around 45-50 pb, which is smaller than the present 
calculation. The reason for a lower timelike cross section 
in Ref. fl7j is the higher order contributions that involve 
gluon GPDs were not considered. 

IV. SUMMARY 

We have investigated the possibility for the inclusive 
and exclusive dilepton detection in the photoproduction 
process at the proposed Deep Inelastic Electron-Nucleon 
Scattering at the LHC (LHeC) machine. The photon- 
proton cross sections have been computed for inclusive 
photoproduction at NLO accuracy and are of order the 
unit of picobarns. The number of events is evaluated to 
be of order 10 5 assuming an integrated luminosity of 10 
fb _1 , which means that such a measurement is fairly fea- 
sible. It was shown that that process is a good place to in- 



vestigate constraints to the parton distribution functions 
for the proton and the photon at very small- x region. 
We have investigated also the exclusive dilepton photo- 
production using the LHeC design. It was found that the 
spectrum is dominated by small invariant masses and the 
forward amplitude scales with (Q 2 at /M 2 ) modulo loga- 
rithmic enhancements, where Q sa t is the phenomenolog- 
ical parton saturation scale. Thus, the differential cross 
section reads as da/dM 2 tx (M 2 )~ 3 [l+log(M 2 )]. The in- 
tegrated cross section gives er CX c(7 +P — > £ + £~ + p) — 80 
pb with the cut M t + t - > 1.5 GeV/c and also one has 
a(ep — > £ + £~ + p) = 2.1 pb for the LHeC energy regime 
and the number of events is around 1.2 x 10 4 . It can be 
considered a lower bound for predictions as the timclikc 
photon configuration would give cross sections bigger by 
a factor three or four. The number of events is one order 
of magnitude smaller than for the inclusive case whereas 
the experimental detection is improved due to the pres- 
ence of one rapidity gap in final state. 
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